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Fractionalized phases are studied in a low energy theory of exciton bose condensate in a muhi- 
band insulator. It is shown that U(l) gauge theory with either fractionalized boson or fermion can 
emerge out of a single model depending on the coupling constants. Both the statistics and spin of 
the fractionalized particles are dynamically determined, satisfying the spin-statistics theorem in the 
continuum limit. We present two mutually consistent descriptions for the fractionalization. In the 
first approach, it is shown that fractionalized degree of freedom emerges from reduced phase space 
constrained by strong interaction and that the U(l) gauge field arises as a collective excitation of 
the low energy modes. In the second approach, complimentary descriptions are provided for the 
fractionalization based on world line picture of the original excitons. The emergent gauge structure 
is identified from the fluctuating web of exciton world lines which, in turn, realizes the string net 
condensation in a space-time picture. 



I. INTRODUCTION 

Fractionalization is a phenomenon where a microscopic 
degree of freedom in many-body system decays into mul- 
tiple modes by strong interaction. For example, it is well 
known that interacting electrons in 1+lD decay into spin 
and charge density waves 0. In 2-1- ID, electrons under 
strong magnetic field decay into quasiparticles carrying 
only fraction of the electron charge^. Given the nov- 
elty of the phenomena, it is natural to ask whether the 
fractionalization can occur more ubiquitously in other 
systems. 

Possible fractionalized phases have been extensively 
studied in 2-f ID quantum spin system in view of its rele- 
vance to high temperature superconductivity (a compre- 
hensive review for this approach can be found in Ref.jj])- 
Earlier, Anderson proposed that a neutral spin-1/2 ex- 
citation, dubbed as spinon, can arise in a quantum dis- 
ordered (spin liquid) phase of 2-l-lD antiferromagnetj^. 
In gauge theory picture various mean-field states have 
been found which have the spinon excitation along with 
a Z2 gauge fieldH S S i- 

Since the Z2 gauge field 
is discrete, the deconfinement phase can be stable with 
massive or massless spinon excitations. If gauge group is 
continuous like U(l), the deconfinement phase with mas- 
sive spinon is not stable at zero temperature in 2+ ID 
owing to the proliferation of instantonsjoj. At zero tem- 
perature, the deconfinement phase can occur if there is a 
time-reversal s ymm etry breaking|lOj. or there are mass- 
less modes [TTllT^ Il3| . The massless modes can occur 
at a quantum critical point or they can occur more 
generically protected by lattice symmetry!^ The 
existences of the deconfinement phase have been also 
demonstrated in an exact soluble spin modelOIII and 
a dimer model in frustrated lattice[3- More recently, 
various 3-1- ID models have been constructed which show 
deconfinement phases[30|. These models have played 
an important role in studying the existence of the frac- 



tionalized phase. However they usually involve spin de- 
grees of freedom defined on bond centers of lattices or 
four spin interaction which break the rotation symmetry 
in spin space. At the same time bosonic models have 
been constructed which may be realized in principle in 
Josephson junction arrays^^i ES- clearly of inter- 

est to search for more examples of realistic models which 
support fractionalization. 

In a sense the fractionalization is opposite to the phe- 
nomena where multiple particles are bound to form a 
composite particle at low energy, such as molecules made 
of atoms and proton (neutron) made of quarks. The lat- 
ter phenomena look more natural in the sense that larger 
structure is constructed from smaller building blocks. 
Fractionalization seems to work in the opposite direction. 
The puzzling nature of it becomes sharper when we real- 
ize that the fractionalization is about low energy physics. 
Fractionalization concerns the emergence of atom-like ob- 
jects out of 'molecules' at low energy, not probing real 
'atoms' inside 'molecules' at high energy. This concep- 
tual difficulty of identifying fractionalized degree of free- 
doms out of original particle translates into the difficulty 
in the formalism to find a unique way of decomposition. 
There is no unique way of assigning statistics and spin to 
fractionalized particles. For example, spin liquid phase 
can be described either by bosonic spinon's''! or fermionic 
spinon i^J. Why should we choose one way of decom- 
position over another ? More importantly, what is the 
fractionalized particle (emergent atom) after all, if it is 
not the real constituent (real atom) inside the original 
particle (molecules) ? 

The purpose of the present paper is to answer these 
questions by demonstrating that not only the emergence 
of fractionalized particle but also the statistics and spin 
of the fractionalized particle are all consequences of the 
intrinsic dynamics of the model. We explicitly show that 
both fractionalized boson and fermion can emerge out 
of a single theory depending on coupling constants, and 
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that the spin and charge of the fractionaHzed particles 
are also uniquely determined from dynamical constraint. 
The microscopic system we consider is a model describing 
excitons in multi-band insulator. 

Since fractionalization is a rather unfamiliar concept 
in condensed matter physics, the question is often asked, 
how is it possible to visualize the new particles and the 
gauge field in terms of the original particle ? In this paper 
we provide an answer by constructing a world line rep- 
resentation of the fractionalized particles and the gauge 
field in terms of the world lines of the original particles, 
which are in our case excitons. The world line (or dual- 
ity) picture is complimentary to the more common field 
theoretical derivation, in that the confined state which 
is difficult to treat in field theory language, is easier to 
visualize. Furthermore, a qualitative understanding of 
the transition between confined and deconfined phases 
is possible in terms of pictures. It is our hope that by 
explicitly working out this example using both methods, 
we can gain deeper insights into these fascinating new 
phenomena. 

We consider a system of spinless fermions with N de- 
generate bands where each band is composed of a con- 
duction band and a valence band. The Hamiltonian is 

N 

N 
a—l k 

a,f> k,k' ,q 

Here (d^) is the annihilation operator of the fermion 
in the a-th conduction (valence) band, e^. (e^), the en- 
ergy dispersion of the conduction (valence) band, /i is 
the chemical potential. For simplicity we set — ii — 

A= 2^ and - A* = -^"^ ~ The band gap is 

^ ^ jg tuned so that there is equal density of 

conduction particle and valence hole. The index a labels 
the band degeneracy which is N for both conduction and 
valence bands. We imagine that the index labels orbital 
degeneracy, but we shall refer to it as the flavor index 
following the particle physics literature. U is the over- 
lap integral which admixes the conduction and valence 
bands. V^^ is the strength of the Coulomb repulsion 
between fermions in band a and band h. It is assumed 
that the Coulomb interaction is local in real space and 
that the strength of the interaction depends on whether 
fermions are in a same band or in different bands, that 
is, y^"'' ^ V + V'Sab- V = L^-i is the volume of the 
system in d space dimension d — {D ~ 1). We single 
out the attractive interaction between conduction parti- 
cle and valence hole. The rest of the Coulomb interaction 
is assumed to be already included in the band structure 
determination of A'^ + A'^ by the Hartree-Fock approxi- 
mation. 



The attractive interaction leads to a bound state be- 
tween particle and hole which is an exciton. If the inter- 
action energy is smaller than band gap, the exciton state 
lies between the gap. Since it is metastable, the fermion 
and hole will eventually rccombine in a finite time. If the 
binding energy is large enough to overcome the band gap, 
particles and holes are spontaneously created to form a 
bose condensate in the ground state. We will consider 
the latter case. 

To fix ideas, it is useful to consider one particular real- 
ization of this model. Consider a hetero-junction where 
two thin semiconductor layers are separated by a thin 
barrier, shown schematically in Fig. ^ The energy 
_|_ |-]-^g separation between the bottom of the 

conduction on the left hand side and the top of the va- 
lence band on the right, and can in principle be tuned by 
applying a voltage. U is the tunneling matrix element of 
an electron across the barrier. In our model only elec- 
trons in bands with the same orbital index are allowed 
to tunnel, which is reasonable if we assume that different 
orbitals are orthogonal to each other due to their symme- 
try within the plane. If the Coulomb attraction between 
an electron on the left and a hole on the right is strong, a 
finite density of excitons will form spontaneously in the 
ground state, resulting in a Bose condensate of excitons. 
In real experiments it has proven difficult to realize the 
Bose condensation of excitons in the bulk, because the 
excitons tend to have attractive interaction, leading to 
the formation of exciton molecules or electron-hole liq- 
uid ground states. The most promising current example 
is indeed the layer geometry shown in Fig. in GaAs 
hetero-junctions. Instead of using a bias voltage, the 
conduction and valence bands are populated by optical 
pumping, leading to a meta-stable state|23. The dipo- 
lar interaction between the excitons assures that they are 
repulsive, making the system a promising candidate for 
exciton condensation. Here we are mostly interested in 
Eq. Q as a model system which we assume can be re- 
alized in both two and three dimensions. The presence 
of N degenerate bands lead to additional complication in 
that there can be iV^ species of excitons and some or all 
of them may or may not condense. The order parameter 
for exciton condensation is 

r" c? >, (2) 

where i is a site index. In this paper we examine the 
case where = Xd 0, i-e., the diagonal exciton are 
condensed. While it is possible that the symmetry break- 
ing happens spontaneously, the presence of the U term 
in Eq. explicitly breaks the symmetry and assures 
that ^ 0. (In the case where the condensation of di- 
agonal excitons occurs spontaneously, we can set U = 
in our model. ) We also assume that there is a finite 
amplitude for off diagonal exciton condensation. In ap- 
pendix A we work out the mean field theory to show 
that this assumption is locally stable in some parameter 
range. In the remainder of the paper we examine the ef- 
fect of phase fluctuations of x"** = Xos'^ , to see whether 
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FIG. 1: Schematic diagram of band structure in the hetero- 
junction of semiconductors. 



the ofF-diagonal condensate may be destroyed, and to see 
under what circumstances a novel state of matter with 
fractionahzed particle and gauge field may emerge. 

Fig. 121 provides a flow chart for the development of 
this theory and can be used as a road-map to this paper. 
Starting from the microscopic Hamiltonian given by Eq. 

we can derive a Landau theory for the order parame- 
ter x"'' based on symmetry considerations. This is done 
in section II and the resulting Landau theory is shown in 
Eq. l(TH) . Since x"" are assumed to be condensed, the re- 
maining degrees of freedom are just x""^ — Xoe'^° and we 
focus only on phase fluctuations. Note that x""^ and x^" 
are in general independent because x""^ is the pairing of 
conduction electron in band a with valence hole in band 
6, which is distinguishable from pairing electron in band 
b and hole in band a. Nevertheless, due to the condensa- 
tion of x°°, there is a term which locks the phases, which 
is proportional to ^2 Ea#b cos(6lf + 9^"). For K2 < 0, 
0f = -e\'' is favored. Thus x"^ is locked to x*""* and 
the Landau theory reduces to a Hermitian matrix model 
Eq. ini). Similarly, for ^2 > Of = -O^^ + ttis favored 
and we have the anti-Hermitian matrix model shown in 
Eq. (|18|l . Further development depends on the next or- 
der term in Xo- For the Hermitian model, the third order 



term K'^ cos( 
constraint 



- 9^ -I- O'l"-) is allowed which leads to the 



= 



(3) 



for strong coupling. This constraint can be satisfied by 
writing 6'f'' = 4)°^ — (t)\. The field e"^> emerges as the bo- 
son field which carries the flavor quantum number, but 
it is coupled to a U(l) gauge field. In d = 3 and for 
large enough iV, a deconfined phase is possible. It is 
called the Coulomb phase, because it mimics our world 
of non-compact QED coupled to charged particles. The 
Coulomb phase lies between the confined phase and the 
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FIG. 2: Flow chart from the microscopic exciton model to 
low energy theories showing various fractionahzed phases. 



Higgs phase as a function of coupling constant k (see Fig. 
131(a)). The Higgs phase is the phase where the off diago- 
nal excitons are condensed and the confined phase is one 
where the excitons are disordered by phase fluctuations, 
with exponentially decaying correlation functions. The 
Coulomb phase is the novel phase which exhibits gapped 
fractionahzed bosons and gapless U(l) gauge field which 
are analogs of photons. This is discussed in section HI. A. 

In section V. A we discuss this model again using the 
world-line representation. It is natural that the off di- 
agonal exciton x"'' is represented by a double line with 
a, b labels. In the strong coupling limit, these excitons 
scatter with each other and exchange partners rapidly 
(see Figs. 01and|Sl). Thus the individual exciton loses its 
identity on a short scale ^1. However, a single line car- 
rying a flavor index may survive over a long scale as 
shown in Fig. (HI These single lines represent the world 
lines of the fractionahzed particles. However, it is clear 
that these world lines are not freely propagating, but 
connected to each other via a web representing multiple 
scattering with other excitons. The web of excitons turn 
out to show the same dynamics as the world sheet swept 
out by a line of electric flux in a compact U(l) gauge the- 
ory. This is reviewed in section IV. The size of the web 
gives rise to a third length scale ^3 shown in Fig. O The 
divergence of ^3 signals the transition between the con- 
fined phase and the Coulomb (deconfined) phase. Finally 
the divergence of ^2 leads to the Higgs phase of the cou- 
pled boson-gauge field system, which corresponds to the 
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Bose condensation of the original off diagonal excitons. 
Next we turn to the anti-Hermitian model, where 
= by symmetry. The next most important term 
is K4COs{9f + e^" + ef + ef'^). For ^4 < the term 
favors 



nhc 



Sf" = 0. 



(4) 



These constraints are less restrictive than Eq. ©. How- 
ever the set of equations I0J are not independent and it 
can be shown that they are satisfied by 



Qa6 



(5) 



where Zi = ±1 is an Ising degree of freedom on each site. 
The resulting theory consists of the boson e'*^' coupled to 
a U(l) gauge field plus an Ising degree of freedom which 
is decoupled at low energies, as shown in section III.B. 
The world line version is discussed in Sec. V. B. 
Finally for K/^ > 0, the constraint 



Qab 



Qbc 



Qcd 



nda 



(6) 



can not be satisfied by any bosonic parameterization. 
This leads to a frustration among the phases which sup- 
presses the condensation of off diagonal excitons. How- 
ever, a fermionic decomposition satisfy the constraint 
perfectly. A U(l) gauge field is also generated, but it 
turns out that a background 7r-flux is the favored saddle 
point. This leads to the familiar Dirac fermions centered 
at (±^, in 2-l-lD. The low energy fermion degrees of 
freedom is represented by a four component Dirac field, 
which can be interpreted as carrying spin degree of free- 
dom, in agreement with the spin-statistics theorem. This 
is detailed in section III. C. It turns out that the fermion 
also emerges from the world line picture, because we can 
show that each closed loop of single world line is asso- 
ciated with a minus sign, as is shown in section V. C. 
However, the spinor nature of the fermion requires a de- 
scription of the short distance physics which is not easy 
to capture in the world line picture. 

The readers who are interested only in the low energy 
physics of the phenomenological model may skip the first 
part of Sec. II and start from the Landau model, (|14|l . 
The readers who are familiar with the dual representa- 
tion of the compact U(l) gauge theory may skip the Sec. 
IV. Sec. HI and Sec. V deal with the same phenomena in 
two complimentary approaches . Sec. HI provides a more 
quantitative but abstract analysis of the phenomenon of 
fractionalization. Sec. V is less quantitative but more 
physical insights can be obtained on what is really oc- 
curring in the original exciton language. 



II. EFFECTIVE THEORY OF EXCITON BOSE 
CONDENSATE 

Starting from Eq. we treat the attraction be- 

tween the conduction electron and the valence hole by a 



Hubbard- Stratonovich decoupling. The effective action 
for the exciton condensate and the fermions are obtained 
to be (for details, see Appendix A), 



S = S1 + S2, 



(7) 



where 



k a 

q a.b \ ^ k 

~^T.dt{x\rci^X (8) 



and 



^2 = dT 



VV 

v + v 



7 EE 141' 



q a 



(9) 



Here Xk is the N x N complex matrix for the exciton con- 
densate field in momentum space. Each component x'k' 
describes the condensate of excitons made of particle in 
the a-th band and hole in the &-th band. (Note that or- 
der of indices are important, x"'' 7^ X^°'-) Self-consistent 
equations are solved for x"**- Various solutions are pos- 
sible depending on the values of the V, V and U. In 
Appendix A, a few saddle point solutions are obtained 
for = 3. In the following we concentrate on one of the 
solutions, 



at 



X 



Xo + iXd - Xo)Sab, 



(10) 



where both the diagonal and off-diagonal components 
have finite amplitude. 

With the exciton condensation the single particle ex- 
citation gap becomes larger than the bare band gap be- 
cause it takes extra energy in breaking the particle-hole 
binding. At energy far below the single particle excita- 
tion gap the fluctuations of the exciton bose condensates 
become effective degrees of freedom. Integrating out the 
fermions we obtain the effective action for the exciton 
condensate fields x- Symmetry enables us to construct 
the low energy theory without actual calculation. Since 
the fermions appear only in Si the symmetry of 6*1 de- 
termines the form of the action generated by their in- 
tegration. Si is invariant under global U{N) x U{N) 
transformation. 



Cfc Uck, 
Xk ^ UxkV\ 



(11) 
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where 



(dldl. 



,0 are N- 



component vectors for the fermion fields and U and V are 
independent U (N) matrices. The symmetry put strin- 
gent constraint on possible term in the effective action 
because it must be invariant under the last line in Eq. 
()f f |l . Up to the sixth order we write the resulting effec- 
tive action, 

51 = ^ y dr krtrxldrXi ~k ^ tr{xl_^_fXi + h.c.) 

i f=i,...,d 

+k2trx\xi + k^{trx\xif + kitrx\xix\xi 
+kQtrx\xix\xix\xi + ^6(*''XiXi)(*''XiXiXiXi) 
+h{trx\xif\ (12) 

Here Xi is the exciton condensate field at space coordi- 
nate i. f is a unit vector in the d-dimensional square 
(cubic) lattice for d = 2 (d = 3). In principle, all of 
the parameters k^., k and A:„ can be determined from the 
Hamiltonian Note that there is only linear deriva- 
tive in the imaginary time direction because the exciton 
is non-relativistic boson. It is convenient to discretize the 
imaginary time direction to rewrite the above action as 



-S'l = XI [ " yi'^xI+ArX* - X tr{x\+.x^ + c.c.) 

i r 

+K2trxlxt + K^{trx\xif + Kitrx\xix\Xi 
+K6trxlxixlxixlxi + Kf^{trxlxi)itrxlxixlxi) 
+K';itrxlx^f]■ (13) 

Here i — (r, i) is index for lattice point in D-dimcnsional 
Euclidean space-time. The dimensionless coupling con- 
stants in the lattice action are given hy K = 2kr, 
K2 = ATfc2 + kr and Kn — Arfc„, where the discrete 
time step is set to be Ar = || . 

The term ^2 in Eq. Q has lower symmetry. The first 
term is proportional to V which distinguishes between 
the interaction between particles with the same and dif- 
ferent flavor indices. The U{N) x U{N) symmetry is 
broken down to C/(l)^^ where each U(l) is generated by 
diagonal U{N) matrices U and V, i.e., Uab = (^ahc"^" and 
Vob = SabS^'^ ■ The second term is proportional to U 
which explicitly break the phase symmetry of the diago- 
nal exciton condensate. This further breaks the symme- 
try down to [/(l)^ where each U{1) is generated by the 
diagonal matrix hi — V . In discrete space time the total 
action is easily written as 



S 
S2 



Si + 5*2, 



(14) 



El 



+xr) 



(15) 



where Ky> = ^ 



-At and Kn — 



vu 



v+v 



tAt. 



The U{1)^ transformation rotates only the phase of 
each component x"''. ( Strictly speaking, the action l|14|) 
has only U{1)^^^ symmetry because x itself is invari- 
ant under one U(l) transformation with U — V = e**^.) 
Thus all amplitude fluctuations of x"** are gapped and we 
ignore the amplitude fluctuations. With the nonzero am- 
plitudes for all exciton flelds at the saddle point I|1U|) . we 
have N"^ independent U(l) phase modes. The phases of 
the N diagonal modes x"'"' are fixed by the explicit sym- 
metry breaking term Kjj. Thus the fluctuations of the 
diagonal modes are gapped. Alternatively, in the case 
where the symmetry breaking of the diagonal excitons 
is spontaneous and U — 0, there will be N Goldstone 
modes. In the following we will focus on the dynamics of 
the remaining N'^ — N off-diagonal modes which will be 
the same in both cases. 

The effective action for the off-diagonal modes can be 
readily obtained from the full action p4|l . The phase 
locking of the diagonal modes plays two important roles 
here. First, it create mass for some of the off-diagonal 
modes. We pick out terms in Eq. (|14|l which are propor- 
tional to Xo and products of the the diagonal elements 
X°°. The derivation is straightforward and only the re- 
sults are shown here (for details see Appendix B). This 
is analogous to the Higgs mechanism. Depending on the 
sign of the 'mass' term for the off-diagonal modes, 



k2xij:^ 



nab 



nba \ 



(16) 



with K2 — 2x 



K, + NK', {xl + {N~ l)xl) + SKexl 
X is constrained to be either Hcrmitian or anti-Hermitian 
at low energy. Second, the coherent diagonal modes make 
the dynamics of the off diagonal mode relativistic. Once 
X is constrained to be Hermitian or anti-Hermitian ma- 
trix, the excitons x°^ and x''" are no longer independent, 
but they are anti-particles to each other. The physi- 
cal picture is that two excitons x"*" and x^" scatter with 
each other to become two diagonal excitons x"" and x''^- 
The diagonal excitons are condensed into vacuum and 
from the point of view of the off diagonal excitons, this 
is an annihilation process of two excitons. Conversely, 
two off diagonal excitons can be created out of vacuum. 
Since the low energy theory is relativistic, the effective 
action has symmetric form in space and time directions. 
It reduces to a Hermitian matrix model (for details see 
Appendix B), 



5 = - J E tr'{xlx,+ h.c.) + KnY.tr' xl 

n>3 i 



(17) 



for K2 < or an anti-Hermitian matrix model, 
S = -^Y. i'^'(xlxj+/i-c.) + E^2„Zir'; 

<i,j'> n>2 i 



(18) 
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for K2 > 0. Here i,j are indices for lattice in D- 
dimensional Euclidean space-time, x is Hermitian (anti- 
Hermitian) matrix without diagonal element in 117() 
(dl))- The ' sign in the trace is to recall that there 
is no diagonal element in x- The explicit form of Kn is 
also obtained in the Appendix B. The first term in the 
above actions is the kinetic energy term for the cxciton 
condensate and the second term, the potential energy 
term. 

Some remarks are in order for these models. First, 
(|17|l and H18|) are two different theories in general. How- 
ever they are equivalent in some special cases. The Her- 
mitian matrix model with K2n-i = and = 
is equivalent to the anti-Hermitian matrix model with 
K2n = (— l)"fc2n. They are related to each other by 
the transformation x — *X- Similar equivalence also 
exist within the Hermitian matrix model. The Hermi- 
tian matrix model with Kn is equivalent to the Hermi- 
tian matrix model with = (— 1)"X„. They are re- 
lated by the transformation x = ~X- Second, the lowest 
order nonzero interaction is most important in the sec- 
ond terms of 1)17(1 and p8|l if Xo << 1- This is because 
the effective phase stiffness for the off-diagonal modes is 
suppressed by factor of x". More precisely, nonzero Ki 
with smallest I determines the dynamics of the theories 
if \Ki\ >> x^IKj+toI for all m. The above consideration 
greatly reduces the number of theories to consider. Up 
to the quartic order, we have only the following three 
different low energy theories, 

A. Hermitian matrix model with 7^ 

B. anti-Hermitian matrix model with K^^ < (equiv- 
alently, Hermitian matrix model with ii'a ~ and 

< 0) 

C. anti-Hermitian matrix model with K/^ > (equiv- 
alently, Hermitian matrix model with — and 
K4 > 0) 

(19) 

In this paper we examine the low energy physics of the 
above three cases in detail. We note that the cases in 
parenthesis in cases B and C are non-generic and require 
some tuning parameter to set K3 = 0. We do not con- 
sider non-generic cases with = — where more 
low energy degrees of freedom may emerge. 

III. EMERGENT U(l) GAUGE THEORY IN 
THE MATRIX MODELS 

In this section we will analyze the low energy physics 
of the three theories H19|l in the strong coupling limit 
l^nlXo >> 1- The degrees of freedom of the (anti) Her- 
mitian matrix will be further reduced by the dynamical 
constraints imposed by the strong interaction. We will 
see that the low energy degrees of freedom are U(l) gauge 



field and particles which carry only fractional quantum 
number of exciton. The statistics and spin of the frac- 
tionalized particle will be determined by the coupling 
constants. 



A. Emergence of fractionalized boson and U(l) 
gauge field 

Here we consider the first case in ((19(1 : the Hermitian 
matrix model ((17(1 with if 3 ^ 0. Without loss of gener- 
ality we assume K3 < 0. For simplicity we consider the 
case with Kn < for all n even though only the K3 is 
important for small Xo- The K3 term in l(17() . 

^3Xo + + ^"') (20) 

a,b,c 

leads to constraints (in addition to the Hermitivity) , 

in the strong coupling limit ji^alxo >> 1- These condi- 
tions also constrain sum of n > 4 phases. For example, 
by summing the two conditions 

Qac^Qcd^Qda ^ Q (-22) 

and using the Hermitivity 6*^° + 9°''^ = 0, we obtain 

gab ^ gbc ^ 0cd _^ Qda ^ Q (-23) 

Thus the higher order interaction terms in 1(17(1 are also 
minimized with Kn < 0. The number of independent 
constraints imposed by is (^-iK-'^"^) leading 
to _ iN-i)iN~2) ^ ^ _ 1 goft j^odes. These 

modes can be parametrized by the N phases 0" so that 

0f = ~ 4>\ (24) 

with a redundancy (/)° ~ (/)^ -I- (pi. The (N-1) low energy 
bosons are the soft modes associated with the ?7(1)^~^ 
symmetry of the Landau theory 1(14(1 . Note that e*"^' 
carries only one flavor quantum number while exciton 
carries two. Thus we refer to this boson as slave boson. 
At this stage, the slave boson is just a book-keeping tool. 
It may or may not appear as low energy excitation. If 
the slave boson appears as true low energy excitation, we 
will refer to it as fractionalized boson, the term which is 
mainly used in literatures. 

For the low energy modes 0° there is no potential en- 
ergy since Eq. ((20(1 is minimized. The effective action 
becomes 

^ = "lE E [e'^^^"^^°^e-'(*'-^.') -I- c.c] ,(25) 

a<b <i,j> 
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where k = 2Kxo- Introducing the Hubbard-Stratonovich 
transformation for the quartic term, 



exp ( Je'W-*")e*W?-*?) + c.c. 



xexp(-J[h^''P + h|i'p 

'Hi ^ 'Iji ^ 



(26) 



and parameterizing the two complex Hubbard- 
Stratonovich fields by four real variables, 



1 



(27) 



we obtain the action 
S, = '- 



[2|,7^''pcosh(20 



E 

<i,j> ,a<b 



1 ab I ■ / -\-ab ab\ 



(28) 



Among the auxiliary fields, there is only one massless 
modes associated with the U(l) gauge symmetry, 



+ - 



(29) 



We retain the U(l) gauge field and neglect the fluctu- 
ations of other auxiliary fields. In general the action 
(|28|l is not real but reality is restored at the saddle point 
where a\^°'^ and wf- have saddle points in the imaginary 
directionj2^. In our case, by comparing the first and last 
terms in Eq. (|28|l . we note that the only solution where 
the effective hopping is the same for all species of bosons 
is 



Vati 



(30) 



The massless mode is one where otj is independent of 

fly and the low energy action 



the flavor index, i.e., af^ 
becomes 

a <ij> 



rj — rje 



ii4>i-'l'j-°-ij) 



(31) 



where \fjf^\ = rj. This is the familiar problem of com- 
pact U(l) gauge theory coupled with N bosons. The 



bare gauge coupling is infinite. At low energy the gauge 
coupling is renormalized to be finite. For example, by 
expanding in k and integrating out bosons which hop 
around a plaquette, we produce a gauge coupling associ- 
ated with the Maxwellian term on the unit plaquette in 
the D-dimensional (hyper) cubic lattice, ^ with 
t = j{N — 1)?7. In the small g limit the confinement 
phase occurs. The slave bosons do not appear as excita- 
tion. All excitations are gapped. This is the disordered 
phase of the off diagonal excitons. In the large t limit the 
bosons have phase coherence (Higgs phase) . One of the 
N boson is eaten by the massive U(l) gauge field. There 
are remaining A'^ — 1 massless modes. In terms of original 
exciton, this is the phase where the N"^ — N oS diagonal 
excitons are bose condensed. Due to locking with the 
already correlated diagonal excitons, only A^ — 1 modes 
survive as the Goldstone modes associated with the off 
diagonal exciton condensates. This is the expected re- 
sult from the U{1)^~^ symmetry as discussed in Sec. II. 
These are the possible phases in 2-l-lD. In 2-f-lD the de- 
confinement phase may exist only at the critical point 
between the confining phase and the Higgs phase. In 
higher dimension the Coulomb phase can occur between 
the confinement and Higgs phases. In the Coulomb phase 
the photon emerges as a massless collective excitation of 
excitons. The physical meaning of the photon in terms of 
original exciton will be discussed in the Sec. V. A. In this 
phase the bosons are gapped and interact by exchange of 
photons. The schematic phase diagram in 3-l-lD is shown 
in Fig. 121(a). 

Since the fractionalized boson emerges from the exci- 
ton which is charge neutral to electromagnetic field, it is 
obvious that it is charge neutral as well. The fraction- 
alized boson is structureless and it seems obvious that 
the boson has spin 0. By spin we mean that the boson 
transforms trivially under Lorentz transformation in the 
continuum limit. However we need to take into account 
the coupling to the gauge field in order to determine the 
spin of the boson at low energy. Let us follow only one 
species of boson by integrating out the rest of the bosons. 
For small t we obtain 



S = -t 



,i{<j>i-ij>j-aij) 



<i,j> 

-E^^os^c- 

r; '^9c 



+ c.c. 



(32) 



Here the flavor index for the remaining boson is dropped. 
C denotes loop in space-time traced by the integrated 
bosons, be = §(j o, is the flux associated with the loop. 
9c ^ (jv-i)fi.(c) is the generalized gauge coupling for 
the non-local gauge action with L{C), the length of the 
loop C. In the massive phase for the fractionalized boson 
with t << 1 small loop contributions dominate the gauge 
action. Note that the gauge coupling can be made small 
with large N even for small t. In the Coulomb phase the 



gauge field has the saddle point value 



0. Hence 



the background gauge field for the remaining boson is 
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FIG. 3: Schematic phase diagrams of (a) the Hermitian ma- 
trix model with Jfs 7^ and (b) the anti-Hermitian matrix 
model (equivalently the Hermitian matrix model with Kz — 
) in d = 3. The Ising transition line is not shown in (b). 

smoothly varying in space-time. Thus the boson remains 
to be structureless in the long- wavelength limit. The 
continuum Lagrangian for the emergent boson becomes 

C = Y,i,id,^-a,r + ^Y.f^^ (33) 

where t^^ and g are the renormalized phase stiffness and 
gauge coupling respectively. /^^ = d^j^a^, — di,a^ is the 
field strength tensor for the emergent non-compact gauge 
field. The fi-actionalized boson is in the disordered phase 
due to the proliferation of vortices and emerges as gapped 
excitation which has the relativistic dispersion, 

ojk = \/ w^fc^ -f w? , (34) 

where v = yti/io is the velocity with = t^=o, ti = 
t^^Q and TO, the mass gap. The mass gap vanishes at the 



transition to the Higgs phase. If the background gauge 
field had turned out to be rapidly varied in the lattice 
scale, the emergent particle in the long- wavelength limit 
would have been different from the boson in the lattice 
scale. If this were the case, multiple components which 
correspond to spin quantum numbers could emerge in 
the continuum limit from the doubling of unit cell. This 
turns out to be the case for emergent fermions, as we 
shall see later. 



B. Emergence of fractionalized boson, Ising mode 
and U(l) gauge field 

In this section we consider the second case in (|19|l . 
For simplicity, we consider the Hermitian matrix model 
with < for all n. Equivalently, we may use the 
anti-Hermitian matrix model with K2n — (— l)"~^|^2n|- 
Here we use the Hermitian matrix model because it is 
more convenient to make the coupling constant have 
same sign. The fourth order term in (|18|l . 

Kixl cos(r'' e''"^ + r'^ + 6'''") (35) 

a,h.c.d 

leads to the constraints 

^ab Qbc ^ Qcd ^ Qda ^ Q ^3g) 

for |-ftr4|Xo >> 1- Surprisingly, these conditions also con- 
strain the sum of three exciton phases. By summing the 
three equations, 

ef + e^'^ + ef^ + ei" = o, 

and using the Hermitivity, Of + 6]°- = 0, we obtain 

2^Qcb ^ Qbd ^ Qdc^ = 0. (38) 

Therefore the sum of three phases are fixed up to Z2 
degree of freedom, 

Qcb^gbd^gdc = OrTT. (39) 

The constraint for the four phases H36() is less restrictive 
than the constraint for the three phases (|21|l in the previ- 
ous section and an extra degrees of freedom is needed. 
Thus the low energy degrees of freedom are parameter- 
ized by the N bosonic variables and a Z2 variable, 

ef^<i^'^-<fi, + \{i + zi), (40) 

with Zi = ±1. The Hermitivity is satisfied with the de- 
composition. The K2n term with 71 > 3 is also minimized 
because 

cos(6if''=^ +0^'-''^ + ...ef ""!) = 1. (41) 
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Thus one exciton field is decomposed into three degrees 
of freedom. The two U(l) phases are the same modes 
which appear in the previous section. The Ising mode is 
the flavor independent Z2 degree of freedom. It is noted 
that Zi is an independent degree of freedom because it 
can not be absorbed into (/)° or (^^ . This can be illus- 
trated by the relation Zi = — e'*^^" +^ ''+6'=") -^^j^Qgg value 
is not affected by the change of 0". The Zi describes 
the low energy Z2 degree of freedom for a composite of 
three excitons. Another important consequence of this 
relation is that because Zi is fully determined in terms 
of ^'''"s, there is no gauge symmetry under which the Z2 
variable transforms nontrivially. This is contrary to the 
gauge symmetry of the U(l) variables, (/)^ — > (j>°;+(j)i. The 
structure of local symmetry in the decomposition deter- 
mines how the slave particles are coupled to emergent 
gauge field at low energy. There will be U(l) gauge field 
coupled to the U(l) fields, but no gauge coupling to the 
Ising mode. 

The effective action for the low energy modes becomes 



S 



a<b <i,j> 



Pj)Q-'i'i'P'i-'l>j)giii^i-^i) -(- c.c. 



(42) 



where k = 2KXo- We first decompose the U(l) variables 
from the Z2 variables. 



exp 



xexp(-J[|e-p + |e 



!-ab\2 



(43) 



With the parameterization for the Hubbard-Stratonovich 
fields 



the action becomes 



(44) 



E 

<j,j>,a<6 



[2|e,fpcosh(2a;f/) 



(45) 



The action is not Hermitian. However the Hermitivity is 
restored at the saddle point xf^ = ixf^ and ct,"'' = icf"''' 

'- ''J '■J 'J ''J 



with real x^^ and at^^. At the saddle point we have real 

''J '■J 

action, 



5 = J E N^-e. 



where 



<i,j> ,a<b 
eabl 



(46) 



\^fj\e''i^ "~^''^'''^'~'''^^\ AU the fluctuations of the auxiliary 
field around the saddle point is massive. Thus the U(l) 
phase modes and the Ising mode are decoupled at low 
energy. The saddle point values of the auxiliary fields 
are determined from the self-consistency conditions, 



Sz 



(47) 



The solution with ^z,^<p > satisfies the self-consistency. 
That is, positive gives ^e*^^^'"^^') > and positive 
C^, in turn, gives (^e''^'l't'<l>l)'^i't>''-<t>';)'j > q xhus the 

Z2 variable is ferromagnetically coupled with each other. 
The final result is the same as that of the usual mean- 
field decomposition with AB -^<A>B + A<B>-< 
A >< B >. The merit of the above procedure is that 
one can easily identify the low energy degrees of freedom 
beyond the saddle point aPDroximation|23j . 

The quartic term of the U(l) phase modes are decou- 
pled by the same transformation used in the previous sec- 
tion. Finally, we obtain the compact U(l) gauge theory 
coupled with N bosons and the ferromagnetically cou- 
pled Ising field. 



S 



-^E E [e'(*'-*^°-'^-)+c.c. 

a <i,j> 



(48) 



<«j> 



where t = \{N ~ l)^,r; and = ^N{N ~ l)^^. 

The physics of the U(l) gauge sector is the same as 
the one discussed in the previous section. On the other 
hand, the Ising variable undergoes a second order phase 
transition from the disordered phase to ordered phase as 
tz increases. The low energy dynamics of the Ising mode 
is described by the real Klein-Gordon model. The energy 
dispersion becomes 



(49) 



where vi is the velocity and m/, the mass gap. The mass 
gap vanishes only at the critical point. In the ordered 
phase the composite field of three excitons has long range 
correlation. Since the low energy theories are decoupled 
for the U(l) and the Ising sectors, the Ising transition 
can occur within any phase in the phase diagram of Fig. 
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13 (b). If the Ising transition occurs within the Higgs 
phase, the Higgs phase of the fractionahzed boson with- 
out long range correlation of exciton is possible. In this 
phase pairs of off diagonal excitons are condensed, i.e., 
gi(e»''+e<=<iA ^ p ^j^.^ contrast to case A. 



C. Emergence of fractionalized fermion, Ising 
mode and U(l) gauge field 

In this section we consider the third case in p9|l . Here 
we choose to use the anti-Hermitian matrix model with 
^2n > which is equivalent to the Hermitian matrix 
model with Kin = (— l)"|-^2n|- The potential energy 



is minimized when 



Q-l 0,2,1 



COS 



gaia2 1/30213 



(50) 



6*"!°^ +6*^ + ... + 6*^ = TT. (51) 

These become constraints if KinXo" 1- The con- 
dition (|51|l defines the 'low energy valley' in the phase 
space of 0°^. However there is no bosonic parameteriza- 
tion for the low energy modes because of frustration in 
the interaction. This can be argued in the following ways. 
Firstly, (|51|) is not consistent with the anti-Hermitian ma- 
trix. The constraint Of + 6',^'= + Of + 6*^° = tt leads to 
ef + ef" + e'l" + ef = tt for d = b. However this is incon- 
sistent with the anti-Hermitivity condition 6°i^ -I- 0\'^ — it 
which leads to ef+e\°'+e\''+ef = mode 27r. Secondly, 
the conditions H51|l are frustrated among themselves. 
Summing the two equations Of + 6^" + 9f + 9f°- = tt, 
Oi" + e^" + 9f + e\°' = TT and using the anti-Hermitivity 
gabj^Qba ^ e^i^+ef = 77 we obtain 9f+0f''+ef''+ef'' = 
mod 27r. In other words, the set of constraint given by 
Eq. (|51|l is inherently frustrated and can not be all sat- 
isfied by any set of {0°^} where 6'°'' are c-numbers. 

We are thus motivated to introduce anti- 
commutativity. 



(52) 



Here ip is Grassmann fields. The i is introduced to satisfy 
anti-Hermitivity, e*^» = — (^e^^^ ^ . Because of the anti- 
commutativity, all of the constraints can be satisfied, i.e., 

= -1. (53) 

Here the fermionic functional integral is performed with 
the measure, 









1 Ua<b 




-1 







M 



I e 



(55) 



where the integration in the left hand side of Eq. (|54|l is 
understood to be performed over the low energy valley 
in the energy landscape for the exciton phases. How- 
ever the decomposition (|52|l is not complete yet. This 
is because the fermionic decomposition also predicts 
^g*(Si +Si +Si )^ ^ 0, while it is zero in the original ex- 
citon representation. It should be zero due to the Z2 
symmetry {9f Of + tt) of the action (|THJ|. Thus we 
include the extra Z2 degree of freedom. 



(56) 



with Zi = ±1. The situation is similar to the bosonic 
case in the previous section. 

With the decomposition we have the consistency for 
'-) elementary integrations, 

e'""' e-'^""'^ = \Y^ J^a [d^lj''*dij'']Mil}^*ip'''ip''*ip^ = 1, 

(57) 



and the potential energy is minimized. 



-1 



(58) 



for all {ai,a2, ...,a2ri}- At the same time the product of 
odd number of exciton fields has zero expectation value 
because of the fluctuating Z2 degree of freedom. This is 
also consistent with the property of original field. Thus 
these fermionic variables and the Z2 spins are the low en- 
ergy modes of the system in the strong coupling regime. 
Note that the fermionic decomposition is over-estimating 
the expectation value of potential energy because the ab- 
solute value of the left hand side in Eq. H58() should be 
less than 1. For a quantitative agreement one may have 
to introduce four or higher fermion terms in the measure. 
Here we proceed with the measure (|55f) because the four 
or higher fermion interaction is irrelevant at low energy 
for d>2. 

The partition function is written as 



^ = / En^a(*dC*rfV'?)< 



(59) 



where 

■i a 



(60) 



a<b <i,j> 



with K = 2Kxl- The measure of the integra- 
tion can be cast into the conventional form by the 
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transforniation'24l 



(61) 



Note that this transformation docs not respect the orig- 
inal 'complex conjugacy' of i/'f and ■ However it is 
legitimate because = tpf and 4'2 = i^i* can be re- 
garded as two independent Grassmann variables. The 
partition function in the transformed fields becomes, 



z 

with the action 

a i 



(62) 



+ C.C. 



a<b <i,j> 



(63) 



In the following we will drop ' sign in ipf . The procedure 
of decomposing the Z2 spin and the fermionic variables 
is exactly the same as the bosonic case in the previous 
section. The final low energy action is obtained to be two 
decoupled theories, 



^ = EE^rc-^E E [^'""''^r^] 

a i 
— tz ZiZj. 



+ C.C. 



a <i,j> 



(64) 



The Ising variable can have disordered and ordered 
phases as discussed in Sec. III. B. In the following we 
focus on the fermions coupled with the compact U(l) 
gauge field. Note that this fermion is spinless (structure- 
less). The spin-statistics theorem does not apply here 
because there is no full Lorentz invariance in lattice. It 
is in the long wavelength limit much larger than the lat- 
tice scale where the spin-statistics theorem requires half 
integer spin for the fermion because of the full Lorentz 
invariance. In the following we will see how the fermion 
acquires spin at low energy. Similarly to the bosonic case, 
we focus on one species of fermion by integrating out the 
rest of fermions. We obtain the effective action for the 
gauge field. 



i <i,j> 



(65) 



he = §Q a is the flux associated with the loop C . ^ 
(jv_i)fi.(6) is the generalized gauge coupling for the non- 
local gauge action with L{C), the length of the loop C. 
Note that the sign of the gauge action is opposite to the 



one obtained from the bosonic action H65|) . The difference 
in sign comes from the difference in the statistics of the 
integrated matters. This difference in statistics makes 
the difference in the low energy spin structure, confirming 
the spin-statistics theorem. 

In the small t limit, the unit plaquette term is most 
important in the gauge kinetic energy term. Because of 
the opposite sign in the gauge kinetic energy term, the 
saddle point of is '-Tr-flux' phase where every space- 
time plaquette encloses flux tt. Now we derive the low 
energy theory for the fermion in 2 + ID. Generalization 
to higher dimension is straightforward. For the '7r-flux', 
it is convenient to work in a unit cell which is doubled 
in both the x and y directions, but not the z direction, 
with the gauge choice 



^r+z^r 



= i(-l)'--+'-«. 



(66) 



The above 'tt-Hux' has been used in lattice gauge theory 
to put Dirac fermion in lattice using only one-component 
fermionj23. The difference from the Ref. is that in 
our case the z direction is the imaginary time direction. 
More importantly in our case the V-flux' is the dynamical 
consequence while it is assumed in Ref. psf . In the 'tt- 
flux' phase the action for the remaining fermion becomes 



S = 



E E ^*IA^IA 



I A=l 



I 



] [i'4'*i2'<Pn - 
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(67) 



where / is the unit cell index and A = 1,2,3,4 la- 
bels the sites within a unit cell. Introducing four- 
component 'spinor', = (V'n, V|2: V'J3>/4)> *T = 
(-0/1, V'/2, "0/31 ■0/4)1 the action is written 



s = E*^ 



i-t 



-2iA, 







Ay 




2iA, 










A, 


~2iA, 


-iAx 


-A, 







2iA, 



+0(A2), 

where / is identity matrix and Ap^*/ — ^fj 
rescale the coordinates A^, ^ dx, A, 
the continuum limit. Introducing a new basis 



(68) 



dy, 2 A, 



with 



U = 



2t 




1 


0-1 



1 


-1 



, We 
dz in 

(69) 



(70) 
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and allowing the fluctuation of the gauge field, we obtain 
the canonical Dirac Lagrangian, 



£ = ^' 



where 



7m 



cr 



+ 4/'- (71) 



(72) 



with the Pauli matrices cr^ for /i = 1,2,3. This is the 4- 
component Dirac fermion with mass j. Note that we 
obtained 4-component Dirac fermion even though the 
minimal Dirac fermion has two component in (2+1)- 
dimensionl2f||. This is because of the parity symmetry in 
our theory. A mass term involving only the 2-component 
Dirac fermion would have broken parity symmetry p^. 

In (2+1) dimension only confinement phase can oc- 
cur for the compact U(l) gauge theory with the massive 
Dirac fermions. However the mass of the fermion de- 
creases as t increases. Thus it is possible that the Dirac 
fermions become massless at a critical point, stabilizing 
the deconfinement phase [l3|. If t increases beyond the 
critical point, the off-diagonal excitons condense. This 
corresponds to pair condensation of the Dirac fermions, 
which generates generalized mass term ip'^* Mabip'' ■ To 
identify massless excitations in the condensed phase, we 
decompose the phase of the off-diagonal exciton into two 
parts as Of = of^"^ + (0^ - 0,^). Here of'^"^ parame- 
terizes the fluctuations which involves the change in the 



sum of phases 6lf + 



and 



the 



fluctuations which does not change the sum. Because 
of the frustration in Eq. (|50|l . there are flat directions 
in both o'^-^^"'^' and 0° classically. Although not shown 

here, quantum fluctuations lift the degeneracy in 0^^''°''' 
because of singularities in the space of classical vacuum. 
This leaves only the [N — 1) massless Goldstone modes. 
The massless spectrum of the exciton condensed phase 
with K/i > is the same as that of the Higgs phase with 
-fi'4 < 0. Thus we expect that the two phases are contin- 
uously connected. It is, of course, a possibility that the 
two phases are separated by a phase transition. 

It is interesting to note that the massless Dirac 
fermions can arise at the quantum critical point of the 
bosonic matrix model. The critical point is described 
by the algebraic spin liquid 28|. It corresponds to the 
deconflned quantum criticality 12J where massless modes 
at critical point stabilizes deconfinement phase. This is 
the first example where fractionalized fermions occur at 
the critical point as opposed to fractionalized bosons p^. 

The frustration in the interaction is closely related to 
the fermionic nature of fractionalized particle. And the 
spin follows the fermionic statistics. In this sense both 
the emergent fermionic statistics and the spin are the 
results of the frustration in the bosonic matrix model. 
In (3+l)-dimension the confinement to deconfinement 
phase transition can occur. For small t, the confin- 
ing phase arises where the fermions are confined by the 



U(l) gauge field. For large t, the massless U(l) gauge 
field (photon) and gapped fermions emerges as low en- 
ergy excitations. The schematic phase diagram in (3+1)- 
dimension is shown in Fig. |21(b). 

It must be emphasized that the fractionalized fermion 
is a different object from the original fermion of the mi- 
croscopic Hamiltonian J^l- The original fermion is spin- 
less and has no internal gauge charge, but it carries elec- 
tromagnetic (EM) charge. The fractionalized fermion has 
spin and charge for the internal gauge field. The only 
quantum number they have in common is the fiavor in- 
dex. Moreover the fractionalized fermion is particle-hole 
symmetric, while the original fermions are in general not 
in that rn^ 7^ m'^ . 



IV. DUALITY MAPPING OF THE COMPACT 
U(l) GAUGE THEORIES WITH MASSIVE 
MATTER FIELDS 

In the next section, we will describe the matrix models 
in dual representation. The goal of the dual description is 
to provide a more vivid picture of fractionalization and 
the emergent gauge field by describing the same phe- 
nomena in terms of the original exciton language. Before 
doing that, in this section we will review the dual repre- 
sentation of the compact U(l) gauge theory coupled with 
matter fields. This will help to identify what plays the 
role of the fractionalized particle and the gauge field in 
the original exciton picture which will be discussed in the 
next section. 



A. bosonic matter 

In this section we discuss the dual representation of the 
compact U(l) gauge theory coupled with bosonic matter 
field. The action (|32|l is rewritten 



.^_i^(cos6c-l). 



(73) 



Here i is lattice site and refers to D positive direction 
of link, t is the phase stiffness of boson and gc, the gen- 
eralized gauge coupling associated with the flux through 
a loop C . Each loop is assigned a unique orientation and 
be — §fj a. The choice of orientation does not matter be- 
cause of the cos(6c'). The summation of C in the second 
term in Eq. (|73|l is over all possible closed loops. If only 
the smallest plaquette is allowed for the gauge kinetic 
term, Eg. H73|l reduces to the well known Abelian-Higgs 
model 29] . Using the Villain approximation, we write the 
partition function. 



Z = J DaJ J2 E 



p,,(i) = — 00 qc= — 00 
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exp(-^ + m) - 0(i) - af,{i) - 27rp^(i)}^ 



'l,fj, 



c 



^9c 



(74) 



Hubbard-Stratonovich transformation for the quadratic 
terms leads to 



j Da j j Dl^ j Dnc ^ 
exp(-i^(/^Wf 

-z^ + /i) - (f){i) - a^ii) - 27rp,,(j)} 



From the identity X]r!= 



= E„=-oo'5(a;-'^) the 



partition function becomes the sum over integer fields, 
Z = jDaJo^ J2 cxp(-i^a^W)2 



(76) 



Integrating out (j) and a we obtain the constraint for the 
current conservation of boson and the flux conservation 
of the electric flux line, 



5^ A^,(^,(j),0 



S{C' ,i,fj.)n , ,0 



(77) 



Here /^(i) is integer defined on link. Nonzero value of it 
at a given line (i,^) represents that the link is a part of 
world lines of boson(s). Negative value represents world 
line of anti-boson. The first constraint in the second line 
of (|77l) implies that the world line of boson should form 
closed loop. This is the result of current conservation of 
boson, nc is an integer valued 'field' assigned to the loop 
C, i.e., it is defined on the space of loops, the set of all 
possible closed loops. The typical size of the loop is given 
by the inverse mass of the integrated matter field. It is 
convenient to think of the loop as a boundary of a sur- 
face. Then nonzero value of nc implies that the surface 
enclosed by C is a part of a world sheet passing through 
the surface. The value with nc > 1 is regarded as mul- 
tiple world sheets and negative value, a world sheet of 
opposite orientation. The sign S{C,i,iJ,) which appears 
in the second constraint in Eq. H77|) is 1 (—1) if the ori- 
entation of C is along (against) the vector ^ associated 



with the (i,fJ.). This constraint implies that such world 
sheet cannot end by itself. In the absence of boson cur- 
rent {lf_i{i) = 0) the end of one world sheet should be 
joined by the end of another world sheet of same orienta- 
tion. This means that they should form oriented closed 
surface in space-time. In the presence of boson current 
{lfj,{i) 0) the world sheet can end on the world line of 
boson. In lattice gauge theory, nc 7^ represents the 
presence of quantized electric flux on that link. Thus the 
bosons serve as sources of the electric flux lines. The 
second constraint in Eq. (|77|l is the flux conservation 
law. Note that the world sheet is not well-defined in the 
lattice scale because of the uncertainty in determining 
a unique surface enclosed by a loop. We can define a 
smoothly varying surface only in the long distance scale 
larger than the typical length of the loops. This is at- 
tributed to the fact that gauge theory emerges only in 
the energy scale lower than the mass of the integrated 
matter field. 



B. fermionic matter 

Here we consider the compact U(l) gauge theory cou- 
pled with fermion in (|65|l . 



^ Va* - E i cos be - 1 



25^ 



(78) 



Here ±/i denote the nearest neighbor link in all 2D di- 
rections and a-^(i + ^i) — —a^{i). Note that the sign 
of the gauge kinetic energy term is opposite to that of 
(|73|l . We introduce an additional phase tt in the flux 
term to apply the usual Villain approximation and the 
Hubbard-Stratonovich transformation. 



-^(-cos bc-1) 



(79) 



For the matter field we expand in t 

Z^p = 



(80) 



where < A >= J Dip* Dip Ae E.'^.''^'. Because of the 
conservation of matter current, only loop of fermion 
world lines contribute to the partition function. 



= E(-i)™ E 

m=0 Ci,...,C„ 



(81) 
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Here m is the number of fermion loops where each loop is 
represented as Cm- There is a distinction of orientation 
for each fermion loop which is determined by the direc- 
tion of fermion current. Unlike the bosonic case there 
is no ambiguity in counting the number of loops. This 
is because two loops cannot intersect owing to the anti- 
commutativity ij^f = 0. L{Cm) is the length of each loop 
and a is the U(l) gauge flux penetrating the loop. 
The alternating sign with the number of fermion loops 
originates from the anticommutativity of the Grassmann 
field. Using ifT^ and and integrating out the gauge 
field we obtain 



Wl 



/„(i)=0,±l «C 



and ±1, e~^/* is replaced with t and there is a factor 
(—1) for each fermion loop and nc- The minus sign for 
each nc originates from the opposite sign in the gauge 
kinetic energy term. 



V. DUAL APPROACH TO THE EMERGENT 
U(l) GAUGE THEORY IN THE MATRIX 
MODELS 



In this section we dualize the three different theories of 
(|19|l in the world line representation of exciton. The frac- 
tionalization, the emergence of photon in the Coulomb 
phase and the statistics of the fractionalized particle will 
be discussed in the exciton language. 



(82) 



Here m(Z^) is the number of fermion loops for a given 
configuration of /^(i). This is similar to the action with 
the bosonic matter in H77|l except that ^^^(i) can have only 



Duality description of fractionalized boson and 
U(l) gauge field 

Consider the effective action for off-diagonal phases of 
the Hermitian matrix with K„ < 0, 



i ^ a<b 



-E?^E E 



n— 3 



* [ai ,a2 , . . . ,a„ 



[ai,a2,...,a„] 



^1=1 



(83) 



Here 9"^^ = —6^"^ is the phase of the Hermitian matrix. We defined a„+i = ai for notational simplicity, k — 2KXo 
and Kn = An\K„\Xo- The summation [ai,a2, ...,a„] is over all ai, a„ modulo cyclic permutation and the reverse 
of the order of indices. No repeated consecutive indices are allowed, that is, ai ^ at+i. /["i-'^a. - ^an] jg ^ symmetry 
factor which compensate the double counting in the cyclic permutation and the reverse of the indices. For example, 
j[i,2,3] _ y [1,2, 3, 4] _ 2^ g^^^ j[i,2,3,i,2,3] _ 1^ g^.^,^ With thc Villain approximation for both the kinetic energy term and 
the interaction term, the partition function becomes 



z = DO J2 n^=3 
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ExE E /'"^"^'•••""'{E^"""^«"2-9[a.,..,a„](0F 



n— 3 



(84) 



I [ai,a2,...,a„ 



1=1 



Introducing two Hubbard-Stratonovich fields If^'^ii) and ul"^ ' ""' (i) and performing the summation over p'^{i) and 
'Z[ai,..,a„] (*)' the partition function becomes sum over integer numbers for the exciton current and interaction vertex, 
, oo 

E exp(-i ^ (/-(.))^-^5: ^(^'"-""'»)^ 

Z^^(z),y[^i---^"l(z) i,fiM<b n^3 I [ai,a2,...,a„] 

/oo n 
i^0exp(-z Y ifm(^'"'i^+f^)-o''i^)}-^T.ll E i^["^-""k*)E^"'"'"'w)' 



i,fi,a<b 



n—3 i [ai,a2,...,an] 



(85) 



where ^[ai - oan] _ i^^^j[ai,a2,...,a2rt\ ^ Here repre- interaction. The free energy cost of unit current in a link 

sents current of ah excitons and yl'^i---'*"] (i), vertex of 
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of exciton with (a) V'^'' = 1, (b) V""" = 1, 



= 2, 



(d) V' 



1 and (e) 



1. An incoming afe-exciton is 



same as an outgoing fea-exciton because x is Hermitian. The 
direction of the exciton world line attached to vertices are 
drawn with the reference of afe-exciton with a < b. 



1' 2 

(d) 



f 2 

(e) 



FIG. 4: Examples of vertices in the world line representation FIG. 5 



Double line representation of the same vertices shown 
- 1, (b) V^'' = 1, (c) 1^123 ^ 2, (d) 
1. The world line of exciton with 
flavor ab is denoted as two lines where one line of flavor a 
moves in the same direction as the exciton world line and the 
other line of flavor b in the opposite direction. 



in Fig. Elwith (a) 

1^1234 ^ 1 ajjd (g) ^1243 



is - and that of unit vertex at a site is ; „."'"„ i . Inteera- 
tion over the exciton phase leads to current conservation 
condition, 

^Af,lf{i)-Y^ (y^<^bas...ar,] _y[baa3...a„]^ ^ 

fj> n as an 

(86) 

The constraint imphes that the vertex ylii''2- -a„] jg 
a source or sink of the currents. Exciton current of each 
flavor is conserved if y[oii2 - o„] _ q ^yi [aia2...a„]. 
jf ^[aia2- -a„] (•j'j ^j^g gj^g j bccomes a source for cur- 
rents and a sink for the currents 
j'^i+f^i^ ...^ ja.ia„ _ Examples of the vertices are 
shown in Fig. Ela)-(e). The vertices describe interaction 
among excitons. One species of exciton is annihilated 
and other exciton is created at the vertex as constituent 
fermion exchange their partners. Even though the exci- 
ton current of each species is not conserved at vertex, the 
flavor current is always conserved. This becomes more 
transparent if we represent exciton world line in double 
lines (see Fig. I^J. 

In the double line representation the link with — 
1 can be viewed as two world lines of a particle of flavor 
a and an anti-particle of flavor h. Each single line will 
be interpreted as the world line of the emergent slave bo- 
son. Since each flavor current is always conserved, the 
single line should form loop. The single line should rep- 
resent the world line of boson not fermion because the 
phase factor associated with the loop of single line is 1 . 
This is consistent with the earlier discussion in Sec. III. 
A. At this stage the boson is just a 'book-keeping' tool 
for the exciton flavor. Whether this boson occurs as a 
low energy excitation or not is a dynamical question. In 



the weak coupling limit (k„ << 1) the energy cost of a 
vertex is very high. Each exciton maintains its identity 
and the boson is conflned within a single exciton for a 
long time I/ku- The low energy physics is described by 
weakly interacting excitons. This is the disordered phase 
of the off diagonal exciton condensate. It is in the strong 
coupling limit (k„ >> 1) where the bosons have chance 
to be liberated. In the strong coupling limit the vertex 
is cheap and exciton completely looses its identity within 
the discretized time (distance) scale. The scattering oc- 
curs so often that the boson is no longer confined within 
a single exciton. The exciton is not a quasiparticle. On 
the other hand the slave boson can propagate for much 
longer time (distance) than exciton as they keep chang- 
ing their partner. Then the slave boson emerges as low 
energy excitation of the system. Some remarks are in or- 
der for the fractionalization. First, it does not mean that 
an exciton is really separated into two slave bosons. In- 
stead a boson is always part of some exciton. The point 
is that a boson is liberated from a particular exciton by 
frequently changing its partner in the strong coupling 
limit. Then the boson (anti-boson) appears to be free 
in the medium of other anti-bosons (bosons). (However 
the boson is not necessarily a quasiparticle because of the 
coupling to an emergent photon as will be discussed be- 
low.) It is emphasized again that the drastic change from 
the exciton quasiparticle picture to the fractionalization 
is due to the strong interaction. Second, the concept of 
fractionalization can be useful in an intermediate energy 
scale even in the case where there is no fractionalization 
in the low energy limit. In the following we explain these 
in more detail. 
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First we consider what configurations contribute to 
partition function . The world hue of the slave bo- 
son forms closed loop and is accompanied with another 
line with opposite current. The simplest configuration 
is two overlapped loops with current of flavor a in one 
direction and b in the other direction. It describes a pro- 
cess where afc-exciton and 6a-exciton pop out of vacuum 
and annihilate. In this sense the ab and 6a-excitons are 
anti-particles to each other even though they are two in- 
dependent excitons. This is a consequence of the conden- 
sation of diagonal excitons. Two diagonal excitons (aa 
and 66) in condensate collide to become two off-diagonal 
excitons {ab and ba) vice versa. The role of the phase 
coherent diagonal excitons is important even though it is 
hidden behind the dynamical stage. Note that we have 
6°^^ — —O'"^ because of the mass term in (|16|) which is 
generated by the phase coherent diagonal excitons. Gen- 
eral configurations consist of a collection of connected 
webs made of boson world lines (world line web). Each 
world line web defines closed surfaces as depicted in Fig. 
Many such closed surfaces may co-exist and interpen- 
etrate each other. It is noted that the closed surfaces 
made of the world line web are oriented. This is because 
the boson current has direction and only two currents 
with opposite direction can share a link. However, can 
we still define the orientation even if there is a crossing 
of fiavor in a vertex ? The answer is yes. One crossing 
diagram is shown in Fig. [T] (a). The crossing vertex in 
Fig. (a) can be obtained from the non-crossing vertex 
in Fig. El(b) through a continuous deformation involving 
rotation as is shown in the figure. The oriented nature 
of the surface has important consequences for the gauge 
group of the emergent gauge field and the statistics of 
the slave particle as will be shown below. Note that 
there are three independent length scales in the world 
line web. The shortest scale is the length of double line 
segment (^i) which corresponds to the life time of a sin- 
gle exciton. The other is the size of the loop of single line 
(^2) which is the life time of the slave boson. Note that 
there is distribution in the size of the loops of single lines 
from the small loop of size ^1 to the large loop whose 
size is much larger than ^1. The smaller loops of size 
^1 are filling the area between the large loops of size ^2 ■ 
Even though there exist large loops of single lines, the life 
time of a single exciton is kept to be ^1. This is how the 
slave boson can emerge as low energy excitations. The 
largest length scale is the size of the bubble made of the 
W6b (^3). In general we have d < < Cs- In the fol- 
lowing the confinement, Coulomb and Higgs phases are 
described in terms of the dynamics of the web. 

In the small k limit, all length scales are of order of 
lattice scale. As k increases ^2 and ^3 increases. How- 
ever ^1 remains almost same for fixed coupling constants 
(k„ >> 1). For ^2 << the closed surface of the world 
line web is well-defined in the length scale a: ^3. Thus 
the low energy theory reduces to the statistical problem 
of fiuctuating surface. The surface can be regarded as 
the world sheet of oriented closed string. The string ten- 




(a) (b) 

FIG. 7: Crossing vertex as a twisted uncrossed vertex. 



sion is determined by the density of exciton and k. In 
the gauge theory picture the string corresponds to elec- 
tric flux line and the string tension, the gauge coupling. 
The gauge group should be U(l). This comes from the 
orientedness of the surface. In contrast, in the Z„ gauge 
theory there are lines where n currents of same direction 
can meet. Then n surfaces can join at these lines. This 
leads to an unoriented surface. If ^3 is finite there is fi- 
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nite tension for the world sheet of electric flux line. ^3 
corresponds the confining scale. How can we probe the 
confining gauge field ? This can be done by examining 
the potential energy between the two slave bosons. 

Consider creating two excitons of flavors ab and be sep- 
arated by distance r. Then ask the probability amplitude 



where i and r is the (D-l)-dimensional space vector and t, 
the imaginary time. The world line of the 6 and d bosons 
will be present only for a finite time scale near t = 
and Ti respectively. After an initial transient time, we 
are left with boson a and anti-boson c connected by the 
world line web as is shown in Fig. |S1 The end point of the 
web is the slave boson and the web is the electric flux line 
connecting them. Thus we effectively measure the poten- 
tial energy between the slave bosons. If ^2 << r << ^3 
two test bosons separated by r are subject to Coulomb 
interaction, V{r) Inr for D = 3 and V{r) ~ —l/r^~^ 
for D > 3. This is because the world line web fluctu- 
ates wildly in this length scale. The linear potential, 

V{r) ^ r, sets in for ^3 << r << ^ where the electric 

flux line is almost straight. For r >> the energy as- 
sociated with the extended electric flux line exceeds the 
mass gap of slave boson. The boson and anti-boson pop 
out to screen the test bosons and create two separate exci- 
tons. Then the potential becomes flat, V{r) ~ constant. 
This is what we observe in the confinement phase of the 
U(l) gauge theory with massive bosons. In the intermedi- 
ate energy scale 1/^3 << E « 1/^2 the effective theory 
is the emergent pure U(l) gauge field with almost gapless 
photon. It is known that pure compact U(l) gauge theory 
is always confined in 2-1-1 dimension. This means that ^3 
is always finite and in the lowest energy scale E << 1/^3 
every excitation is gapped. The Coulomb phase is not a 
stable phase but can exist only as a cross-over. In space 
dimension greater than 2 the real Coulomb phase can 
arise even in the low energy limit. In the following we 
discuss the Coulomb phase. 

As K increases further, the bubbles grow and over- 
lap. The surfaces of two bubble can re-connect to form 
a large and a small bubbles. In this way an infinitely 
large surface can form and the length scale ^3 diverges 
while ^2 remains finite. The world line web condenses 
(but not the individual boson). The low energy exci- 
tation is the collective excitation of the condensed world 
line web. In a given time slice, the intersection of the web 
surface with the time slice forms strings. In the decon- 
fined phase, these strings are infinite in extent. This is 
a space-time picture for the string- net condensation |30l| . 



for the two excitons to end up with excitons of flavors ad 
and dc after r = ri . We follow only the flavor of a and c 
by tracing out flavor b and d and measure 
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I 



ad d c 





a b be 



r 

FIG. 8: Propagation of two excitons separated by distance 
r in the confinement phase. The end points of the web are 
boson and anti-boson, and the web is the electric flux line 
connecting them. ^3 is the length scale associated with the 
transverse fluctuations of the web. 



In the gauge theory language the electric flux lines are 
the strings which are condensed, and the low energy col- 
lective excitation becomes the gapless 'photon'. This is 
the Coulomb (deconflnement) phase. 

What is the massless 'photon' in terms of original exci- 
ton ? To answer the question we flrst construct the Wil- 
son operator. The Wilson operator creates a closed elec- 
tric flux line in the gauge theory picture. Since the world 
sheet of electric flux line corresponds to world line web in 
the exciton picture, the Wilson operator corresponds to 
creation operator of series of excitons. However in order 
to create low energy excitations, the operator needs to 
create excitons in a special order along the loop, as we 
see from the following argument. The world line web cor- 
responding to low energy excitations need to have large 
number of loops. This is because each loop contributes 
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FIG. 9: (a) Series of excitons created by the Wilson operator. 
The solid (open) circle denotes flavor (anti-flavor) of exciton. 
(b) World sheet spanned by the time evolution of the Wilson 
operator. 



FIG. 10: The flavor current caused by the Wilson opera- 
tor, ab indicates an exciton made of particle in the a-th band 
and hole in the b-th band. The Wilson operator induces the 
motion of boson in the counter-clockwise direction and anti- 
boson in the clockwise direction. 



—InN {N, the number of flavor) to the action. From this 
consideration wc identify the Wilson operator, 

w.,,.....„^ E ^^T'^'^T'^-+<>tr\ (88) 

aia2---an 

This operator creates a series of exciton along a contour 
C defined by the sites ii, ...jin- The neighboring excitons 
have same flavor and anti-flavor (see Fig. |51 (a)). As 
the series of excitons propagate in time they can readily 
form a web as is shown in Fig. (b). The correlation 
between flavors enables excitons to efficiently form a web 
maximizing the number of loops. 

What does the Wilson operator do to existing excitons 
? It creates flavor current. For example if we apply 
W ^ ^ior+or+or+<^r} to the states 

> = e^(«f^+»f+«f+«?')|0> (89) 

with |0 > being a vacuum without any off-diagonal exci- 
ton, we have 



W\^2 > = +»r+«i")|0 > . (90) 

Here we used the fact that once ab and be excitons are on 
a site they decay into ac and bb excitons to lower energy 
by condensing the diagonal exciton. Keeping only the 
off-diagonal excitons we obtain the above result. This is 
shown in Fig. ^| The Wilson operator induces flavor 
current in one direction (Fig. IIUI (a)), and anti-flavor 
current (Fig. 1101 (b)) in the other direction, thus inducing 
net flavor current. induces the opposite current and 
the rotational component of the flavor current can be 
identified as 

£jldx^ ^{Wc-Wl), (91) 

where jf is the flavor current. 

With the excitonic expression for the Wilson opera- 
tor we are readily to discuss the physical origin of the 
emergent photon. In the Coulomb phase there exist in- 
flnitely large bubbles of world line web. This causes long 
range correlation between the Wilson loop. Consider the 
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FIG. 11: The connected diagram contributing to the corre- 
lation function between two Wilson operators. 

correlator between the fluctuations of two Wilson loops, 
(5Wl^5Wc), (92) 

where 8Wc = Wc— < Wc > and Ci, C2 are two loops 
separated by r >> ^2 in the imaginary time direction. 
Since Wc\ (W^^) creates (annihilates) world line web, 
the loops Ci and C2 becomes the boundary of the closed 
surface as is shown in Fig. Only the connected dia- 
gram contribute to the correlator because the background 
value is subtracted. In the confining phase with finite ^3, 
the contribution of the connected diagram is exponen- 
tially small for t >> ^3. In the deconfinement phase, 
there are bubbles of arbitrarily large size which connect 
the Ci and C2 ■ This will leads to power law decay in the 
correlator, which signifies the presence of massless mode. 
This is the emergent photon in the Coulomb phase. In 
this way the emergent photon can be understood in terms 
of the exciton language. In the fractionalized boson pic- 
ture, the emergent photon is the long- wavelength fluctu- 
ation of the flavor currents. The flavor current in Fig. ^| 
can be regarded as if the bosons are really liberated from 
the anti-bosons. On the other hand what is really hap- 
pening is the concerted motion of the neutral excitons 
and the exchange of flavors between excitons on a site. 
Therefore the emergent boson should be electromagneti- 
cally neutral. 

As K keep increasing ^2 grows further. As seen from 
Fig. 13 areas with large ^2 correspond to punctures of 
the bubble. In the length scale a; < ^2 the individual 
loops bordering the punctures on the 2-dimensional sur- 
face begin to be noticeable. The surface is no longer 
closed surface. The individual bosons form the bound- 
ary of world sheets and appear as low energy excitations 




(a) (b) 

FIG. 12: A large loop surrounded by small loops in the ex- 
citon world line web (a). In the length scale comparable to 
the large loop, the small loops are regarded as smooth surface 
connecting the large loop (b). 

(not necessarily quasiparticle). For a given length scale 
X, the loops with sizes smaller than x form smooth sur- 
face and are regarded as parts of electric flux line of the 
U(l) gauge field. The loops with sizes larger than x are 
regarded as the world lines of the fractionalized boson. 
This is displayed in Fig. ^1 The size of the large loop 
represents the inverse mass of the fractionalized boson. 
A finite loop implies a nonzero mass of the fractionalized 
boson. Note that only one flavor of exciton is observable 
in the large length scale because the smooth surface does 
not carry net flavor. This is how the fractionalized boson 
arises. In a given time slice, the boson and anti-boson 
are the end of strings formed by the intersection of the 
web surface and the time slice ^0^. So the low energy the- 
ory becomes the non-compact U(l) gauge theory coupled 
with gapped bosons. In the large k limit not only ^3 but 
also 1^2 diverges. The slave boson itself condenses. This 
is the Higgs phase. Unlike the Coulomb phase the Higgs 
phase can occur even in the 2+1 dimension. 

B. Duality description of fractionalized boson, 
Ising mode and U(l) gauge field 

The duality mapping for the Hermitian matrix model 
with K2n-i — and < is exactly same as the 
one in the previous section except for the fact that only 
vertices of even order are allowed. Thus we will not re- 
peat the same procedure. Instead we will discuss the new 
feature which arises from the absence of vertex of odd or- 
der. As discussed in the Sec. III. B, there exists extra 
Z2 degree of freedom which undergoes the disorder to 
order phase transition. How can we interpret the phase 
transition in the dual representation? We will show that 
the phase transition corresponds to a spontaneous gen- 
eration of the odd order vertex. To see this, we first note 
that the presence of quartic vertex implies the presence 
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FIG. 13: The generation of the sixth order vertex from the 
quartic order vertices. 




FIG. 14: Two sixth order vertices connected by the world 
lines of three excitons. 



of all higher vertices of even order. For example, Fig. 
1131 shows the generation of a sixth order vertex from the 
one-loop contribution. In the length scale larger than 
the loop, the sixth order vertex can be regarded as a lo- 
cal interaction. Two sixth order vertices can effectively 
act as two third order vertices connected by world lines 
of three excitons as is shown in Fig. ^1 In the disor- 
dered phase of the Ising variable, the composite of the 
three excitons has only short range correlation. Thus the 
two third order vertices in Fig. Elis linearly confined to 
each other. In the ordered phase the world line of the 
three excitons can be infinitely large. Thus the third or- 
der vertices are unbound. One can regard the term 
as an explicit symmetry breaking term for the Z2 sym- 
metry of the exciton Lagrangian. Then the generation 
of the third order vertex without bare is the spon- 
taneous symmetry breaking of the Z2 symmetry. In the 
world line picture it corresponds to the binding unbind- 
ing transition of the effective third order vertices made 
of the even order vertices. 
C. Duality description of fractionalized fermion, 
Ising mode and U(l) gauge field 



Here we examine the anti-Hermitian matrix model 
with > 0. The duality mapping can be performed in 
a parallel way to the mapping for the fractionalized bo- 
son except for the two differences. First, the coefficients 
of the potential energy term has opposite sign. Second, 
the matrix is anti-Hcrmitian. These two differences cause 
the transmutation of the statistics of the slave particle. 
The effective action for the off-diagonal phases is written 



i ^ a<b 



^^^n^ ^ j[oi,02,...,02„] 



I [ai,a2,...,a2„ 



COS (E^'°'°''^'(*) - tt I - 1 



(93) 



where 0"'' = —9^°" -I-tt is the phase of the anti-Hermitian matrix. Other notations are same as those in Sec. V. A. Note 
that that there is an extra phase tt in the potential term accounting for the opposite sign of the coupling constant. 

From the Villain approximation and the Hubbard-Stratonovich transformations we obtain the similar partition 
function to H85|) . 
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The extra phase factor Qv m the partition function is the 
consequence of the opposite sign in the coupling constant. 



It is given by 
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2 [ai ,a2 ,. . ■■,0.271 



21 



and it can be simply written as 

Wv — e , 



(96) 



where V is the total number of vertices in the polygon 
made of the exciton world lines. Besides Qv, there is 
another phase arising from the anti-Hermitivity. Since 
9°''' and 9^°' are not independent fields, one has to replace 
6'^" into —9'^^ + TT for a < 5 in integrating over 0°^. After 
integration it leads to the same constraint as H86|l but 
with another phase factor, 

oo 



g[ai...a2„] counts the number of emanating exciton 
world lines attached to the vertex v^°'^'"°'^"-\i) whose 
double index is in descending order (e.g. s^^'^^^^i) = 1 
because there is one descending double index 41 among 
12, 23, 34, 41; similarly we have s^^^^^\i) = 2.) The 
second extra phase factor, in turn, becomes 



(99) 



where 



[ai...a2„]^j-) _ 



^ 2 

1=1 



H+l 



Finally we have 

iab(,)_y[oi...o2„l(i) 



e"(^-^) exp 



(98) 



ere E is the total number of exciton world line seg- 
ments. This is because each exciton world line is attached 
to two vertices and the double index is in ascending order 
in one vertex and in descending order in the other vertex 
(see Fig. [T3)l . 



E (C»)^-EE 



i,fi,a<ib 



f^[ai...a2n] 

ai,a2,...,a2n] 



(y[ai...a2„](j))2^ X 



0,0 



(100) 



The extra phase factor can be related to the number of 
faces F of the polygon because the Euler characteristic 
is a topological invariant, 

XEuler^V-E + F = 2{l-g), (101) 

where g is the genus of the surface. Since g is integer we 
have 

e'""^^-^^ = e*"^. (102) 

On the other hand the number of the faces is same as 
the number of the loops of single line. There is extra 
phase factor (—1) for every loop of single line. Thus 
each single line represents the world line of a fermion. 
For example the diagram in Fig. ^] (a) has 4 loops and 
Fig. ^1 (b) has 5 loops of single line. They contribute 
to partition function with the opposite sign. This is the 
signature of the emergent fermion in the exciton partition 
function. As was mentioned before, the orientedness of 
the closed surface is important for the simple relation 
(|102|l because the Euler characteristic 1)1011) can be odd 
integer for non-orientable surface. It is straightforward to 
obtain the same result from the Hermitian matrix model 

withX2„ = {-ir\K2n\- 

The confinement to deconfinement phase transition 
can be understood as the divergence of the scale of the 



world line web as discussed in Sec. V. A. The difference 
is that the bubbles made of exciton world lines contribute 
to the partition function with oscillating sign as the num- 
ber of single line loops increases. The oscillating sign also 
appears in the partition function for the gauge theory 
coupled to fermions as seen from the phase factor e'^"^ 
in Eq. Thus we see that this case corresponds to 

an emergent gauge theory with the opposite sign of the 
kinetic energy term (Eq. JTHJl). 

The phase transition for the Z2 variable can also be 
understood in the same way as discussed in Sec. V. B. 
The difference is that the ordered phase of the Ising vari- 
able is not related to the generation of the usual Kstr 
term in Ijl?!) . It cancels with its complex conjugate term 
for the anti-Hermitian matrix. The spontaneously gen- 
erated interaction in the ordered phase corresponds to a 
flavor dependent interaction, 

E Kf^x'^^x'^x'" (103) 

with Kf" = -iiTf With the flavor dependent cou- 
pling constant, the complex conjugate term no longer 
cancel each other under the anti-Hermitian constraint. 
Eq. (|103|l contributes to (|93|) as vertex of third order. 
Note that the overall sign of Xg"''''^' does not matter be- 
cause the third order vertex should occur in pair because 
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FIG. 15: a&-cxciton is emanating from the vortex Vi and 6a- 
exciton, from the vertex Vj. Thus there is one e~*'^ factor 
associated with the exciton world line either from Vi or Vj. 



(a) (b) 

FIG. 16: Configurations of oxciton world lines involving (a) 
4 loops and (b) 5 loops of single line. 



statistics, in turn, is shown to be responsible for the gen- 
eration of spin in the continuum hmit in agreement with 
the spin-statistics theorem. The fractionahzation is ex- 
plained in two consistent point of views. One is based 
on more conventional approach where the fractionalized 
degrees of freedom are obtained from decomposition of 
original exciton field into multiple modes as required by 
constraints imposed by strong interaction. The other ap- 
proach is based on the world line representation of origi- 
nal exciton. In the latter approach, a more vivid physical 
picture for the fractionahzation is obtained in terms of 
the original exciton language, while more quantitative 
but abstract analysis is given in the first approach. 

We close by making some comments on the observa- 
tional consequences of the emergent photon and fraction- 
alized particles in the Coulomb phase. The photon is the 
only gapless mode and will dominate the specific heat at 
low temperature, as in d = 3 and with cross-over to 
a small gap to d = 2. The fractionahzation of the exciton 
will greatly affect the spectral properties. The exciton is 
not a well defined excitation. The exciton propagator is 
a product of two fractionalized particles in space-time, 
connected by a strongly fluctuating web, i.e. interact- 
ing via the exchange of photons. If the coupling to the 
photon is weak, it can be ignored and the exciton spec- 
tral function is a convolution of two spectral functions of 
fractionalized particles, and will become extremely broad 
with an unusual line-shape. 
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of the constraint between the number of vertices and the 
double line segments, 

nVn = 2E. (104) 

n 

Finally, the massless fermion can arise if the size of the 
loop diverges. 



VI. CONCLUSIONS 



In this paper we studied an exciton model which sup- 
ports various fractionalized phases with either fractional- 
ized boson or fermion along with emergent photon. Phe- 
nomenological matrix model for exciton bose condensate 
is derived from a microscopic theory. The statistics and 
spin of the fractionalized particles arc determined dy- 
namically. The bosonic fractionalized particle is shown 
to occur more generically. The fermionic nature of the 
fractionalized particle is shown to be rooted to frustration 
in the underlying bosonic matrix model. The fermionic 
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Appendix A. Effective action of exciton condensate 
and the saddle-point solutions 

The objective of this appendix is two-folded. The first 
one is to derive the effective action Q-(0 from the mi- 
croscopic Hamiltonian The second one is to present 
the saddle point solutions for the self-consistent equations 
for the exciton condensate field. 

We rearrange the the attractive interaction as 



q ah \ k 



V + V' 



V 



(A5) 



ab k,k' ,q 

ab k,k' ,q 



From now on we omit the prime in the x drop the 
last constant term. 

Now we consider the mean-field solution for the exci- 
ton condensate. We consider a translationally invariant 
saddle point Ansatz, 



k' ' 



(Al) 



ab r~, ah r 



(A6) 



and introduce Hubbard-Stratonovich fields describing the 
exciton condensate to obtain 



for the exciton condensate and solve the self-consistent 
equations for x"**, 



Cah 



v + v 



v + v 



= EE (^")"^iCi'-5^E<v: 



k 

(A7) 



ab q 

tab* 
'^1 



E '^'k^k+q 



(A2) 



fe 



With the anisotropic Coulomb interaction, V — V 

V 



V Sab and a dimcnsionless exciton field, x"'' = f^e 



effective Lagrangian becomes 

VV' 



b\2 



■EEi 



q ab 



V + V' ^ ^ '^9 

q a 



Note that the saddle point equation for has different 
form from Q because of the rescaling and shift in x"''- 

Here we present the saddle point solutions for = 3 
obtained in the 20 x 20 square lattice. For simplicity we 
use the particle-hole symmetric tight-binding dispersion , 
e^ — = — (cfc — m) = —t{cos kx + cosky — 2) + A with t — 1 
and A = 1. Various solutions are possible depending on 
the parameters ofV,V , U. The foUowings are the saddle 
point solutions : 

• 1/ = 10, y' = 5, 



- E E (-5, + 'k + dT {dr + ei - /.) dt 



k a 

V 



" \ " \ " / \ ^ a* ab -lb , \ ^ ib* ab* a 
^^2^2^[2^''k + qXq dk+2_^dkXq C^+q 

^ q ab \ k k / 

+uYT.^^Tdl + dTci). 



(A3) 



k a 



with Xd = 0.68. 



Absorbing the mixing ([/) term into x by a shift of the 
exciton field 



Xf =X:'+^r^5ab5q,^ 



U 



we finally obtain the Lagrangian 

vv' 



(A4) 



q ab 

vu 



b\2 



V + V 



EEixri 



X = 

with xi = 0.4. 
,V^1Q,V'^ -5, U 



= : 




Xd 




Xd 


") 





Xd ) 


= : 




XI 


M 





XI 


XI 


} 


= 0.3 : 





(A8) 



(A9) 



7V^E(:^o""+xr*) 



v + v 

+ J2T.['k*(9r + el-f,)ct 

k a 

+dt*idr + et-fi)dt 



X = 



Xd -Xi X2 
-X2 Xd -Xi 
Xi -X2 Xd 



(AlO) 



with Xd = -0.06, xi = 0.07 and X2 = 0.39. 
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,V=10,V 



-5. U = 0.7 



ab 



Xd Xo Xo 
Xo Xd Xo 

Xo Xo Xd 



(All) 



with Xd = -0.14 and Xo = 0.27. 



The above solutions are locally stable. However it is 
possible that they are not global minima. One may have 
to add further interaction term to stabilize the desired 
saddle point solution. 



Appendix B. Derivation of Hermitian and 
anti-Hermitian matrix models 



In this appendix we derive the Hermitian H17|l and the 
anti-Hermitian matrix model H18|) from the effective ac- 
tion of the exciton condensate (|14|l . This is done by re- 
placing the diagonal element of x by the amplitude of the 
diagonal condensate, Xd in l|14(l and extract the terms for 
the off-diagonal elements of x- In the quartic interaction 



if 4 

XbaXbcXdcXdai 

a,b,c,d 



(Bl) 



there are 5 different possibilities of how the diagonal and 
off-diagonal elements are mixed depending whether 2, 1 
or of the x fields are diagonal. This leads to the action 
for the off diagonal phase modes, 



K4Xdtr' ixX + h.c.) 
+2K4Xdtr ixx^'x + h.c.) 
+K4tr X^'XX'^X, 



(B2) 



where tr {AiA2..An) denotes trace for the product of ma- 
trices which involves only the off-diagonal element of each 
matrix, e.g., tr AB — ^^^i, AabBba- Similarly, the six- 
th order interactions also give rise to potential energy for 
the off-diagonal phase modes, 

Lk, = SKexitr' {xx + h.c.) + Kexltr' ix'' + h.c.) 

+9Kexltr' ixxh + h.c.) + SKexltr' {x\^ + h.c.) 

+^Kexltr x'^X^XX + ^K^x^^r x^xX^X 

+3KeXdtr {xX^XX^'x + h.c.) 

+K(itr XX^XX^XXK 
Lj,, = NK'^ixl + {N~ l)xl)trx^XX^X. (B3) 

6 

L ^1 has the same form as Lka and this renormalizes K4 

into <' =K4 + NK'q{xI + {N -\)xl). Summing Lk„ 
Lxfi and L^,' we have 



{Kixl + iKeX^)tr {xX + h.c.) 



+Kexltr' ix"^ + h.c.) 

+(2<xrf + QKeXd)tr' ixxh + h.c.) 
+3Kexltr'{x''x^ + h.c.) 

+^KeXdtr X^'X^'XX 
+(<' + 6Kexl)tr'x^xX^X 
+3KeXdtr {xX^'xx'^X + h.c.) 
+Ketr XX^'XX^'XX^- 



(B4) 



For Xo << 1 the second order term of x is most im- 
portant. We rewrite the second order term in ljB4|l . 



where 



-^tr ixx + h.c.) = K^tr xt > . cos(r'' + 9"'^), (B5) 



K2 = 2{K^xl + SKexi) 



2x1 



K4 + NK', {xl + iN- l)xl) + m^xl 



(B6) 



If K2 < 0, the energy of (jBSp is minimized by 

6f =. -61^. (B7) 

This gives mass to the ^^^^'~^'> symmetric modes with 
Qab _ Qba ^ rpj^^ remaining modes are ^(^^'^^ antisym- 
metric modes with 9°-^ = -9^°-. For \K2\xl » 1' the 
symmetric modes are frozen and x becomes Hermitian 
matrix (x^ = x)- If ^2 > 0, the energy of ljB5|l is mini- 
mized by 

9f = -61,^" -I- TT (B8) 

and X becomes anti- Hermitian matrix for K2X0 >> 1 
(x^ — — x)- These Hermitivity or the anti-Hermitivity 
condition converts the dynamics of the off diagonal exci- 
ton into relativistic one. To see this we rewrite the first 
term in Eq. 113|l as 



2 
K 



rXt+Xi+ArXi), (B9) 



where we used the Hermitivity or anti-Hermitivity condi- 
tion in the last line. From this, we obtain the Hermitian 
matrix model, 



5 



with 



K 



J E tr'ixUj+h.c.) + J2KnJ2tr'x'l, 

<hj> 



n>3 



K3 = 2[2K';xd + lOKex% 

K4 = kI + IhK^xl 

K5 = QKeXd, 

Kg = Ke 



(BIO) 



(Bll) 



25 



for K2 < 0, and the anti-Hermitian matrix model, 

<i,j> n>2 i 

(B12) 

with 

= < + 3i^6Xd> 
Kq = -Ke (B13) 



for K2 > 0. Here < i,j > is the nearest neighbor in 
the D-dimensional (hyper) cubic lattice, x is Hermitian 
(anti-Hermitian) matrix in Eq. I)B10|I f l)B12p ^. Note that 
all odd order terms vanish in the anti- Hermitian matrix 
model because of the anti-Hermitivity, e.g., 

tr ix^XX^ + XX^'x) = tr (x^'xx'' - X^'xx'') = 0. (B14) 
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